Fractional ac Josephson effect in p- and d-wave superconductors 
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For certain orientations of Josephson junctions between two p^-wave or two d-wave supercon- 
ductors, the subgap Andreev bound states produce a 47r-periodic relation between the Josephson 
current / and the phase difference cj>: I oc sin(0/2). Consequently, the ac Josephson current has 
the fractional frequency eV/h, where V is the dc voltage. In the tunneling limit, the Josephson 
current is proportional to the first power (not square) of the electron tunneling amplitude. Thus, 
the Josephson current between unconventional superconductors is carried by single electrons, rather 
than by Cooper pairs. The fractional ac Josephson effect can be observed experimentally by measur- 
ing frequency spectrum of microwave radiation from the junction. We also study junctions between 
singlet s-wave and triplet p^-wave, as well as between chiral p x ± ip y -wave superconductors. 

PACS numbers: 74.50.+r 74.70.Kn 74.72.-h 74.70.Pq 



I. INTRODUCTION 

In many materials, the symmetry of the supercon- 
ducting order parameter is unconventional, i.e. not s- 
wave. In the high-T c cuprates, it is the singlet d x 2_ y 2- 
wave [1]. There is experimental evidence that, in the 
quasi-one-dimensional (Q1D) organic superconductors 
(TMTSF) 2 X [2], the symmetry is triplet [3], most likely 
the p x -wave [4, 5], with the x axis along the conduct- 
ing chains. Experiments indicate that Sr 2 Ru04 has the 
triplet chiral p x ± ip y -w&ve pairing symmetry [6] . 

The unconventional pairing symmetry typically results 
in formation of subgap Andreev bound states [7] on the 
surfaces of these superconductors [8] . For d- wave cuprate 
superconductors, the midgap Andreev states were pre- 
dicted theoretically in Ref. [9] and observed experimen- 
tally as a zero-bias conductance peak in tunneling be- 
tween normal metals and superconductors (see review 
[10]). For the Q1D organic superconductors, the midgap 
states were theoretically predicted to exist at the edges 
perpendicular to the chains [11, 12]. In the chiral su- 
perconductor Sr 2 Ru04, the subgap surface states are ex- 
pected to have a chiral energy dispersion [13]. Their con- 
tribution to tunneling is more complicated [14] than a 
simple zero-bias conductance peak found for the midgap 
Andreev states. Various ways of observing electron edge 
states experimentally are discussed in Ref. [15]. 

When two unconventional superconductors are joined 
together in a Josephson junction, their Andreev surface 
states hybridize to form Andreev bound states in the 
junction. These states play an important role in the 
Josephson current through the junction [16]. Andreev 
bound states in high-T c junctions were reviewed in Ref. 
[17]. The Josephson effect between two Q1D p^-wave su- 
perconductors was studied in Refs. [18, 19]. Andreev re- 
flection [20] at the interfaces between the A and B phases 
of superfluid 3 He was studied in Ref. [21]. However, An- 
dreev bound states were not mentioned in this paper. 

In the present paper, we predict a new effect for 
Josephson junctions between unconventional nonchiral 



superconductors, which we call the fractional ac Joseph- 
son effect. Suppose both superconductors forming a 
Josephson junction have surface midgap states originally. 
This is the case for the butt-to-butt junction between two 
p x -wscve Q1D superconductors, as shown in Fig. 1(a), and 
for the 45°/45° in-plane junction between two d-wave su- 
perconductors, as shown in Fig. 5(a). (The two angles 
indicate the orientation of the junction line relative to 
the b axes of each d x i_ y i superconductor.) We predict 
that the contribution of the hybridized Andreev bound 
states produces a 47r-periodic relation between the su- 
percurrent I and the superconducting phase difference 
(j): I oc sin(c/>/2) [22]. Consequently, the ac Josephson 
effect has the frequency eV/h, where e is the electron 
charge, V is the applied dc voltage, and h is the Planck 
constant. The predicted frequency is a half of the con- 
ventional Josephson frequency 2eV/h originating from 
the conventional Josephson relation I oc sin </> with the 
period of 2tt. Qualitatively, the predicted effect can be 
interpreted as follows. The Josephson current across the 
two unconventional superconductors is carried by tunnel- 
ing of single electrons (rather than Cooper pairs) between 
the two resonant midgap states. Thus, the Cooper pair 
charge 2e is replaced the single charge e in the expression 
for the Josephson frequency. This interpretation is also 
supported by the finding that, in the tunneling limit, the 
Josephson current is proportional to the first power (not 
square) of the electron tunneling amplitude [23-25] . Pos- 
sibilities for experimental observation of the fractional ac 
Josephson effect are discussed in Sec. V. A summary of 
this work is published in the conference proceedings [26] . 

The predicted current-phase relation I oc sin(c/>/2) is 
quite radical, because every textbook on superconductiv- 
ity says that the Josephson current must be 27r-periodic 
in the superconducting phase difference cA [22]. To our 
knowledge, the only paper that discussed the 47r-periodic 
Josephson effect is Ref. [27] by Kitaev. He considered 
a highly idealized model of spinlcss fermions on a one- 
dimensional (ID) lattice with superconducting pairing 
on the neighboring sites. The pairing potential in this 
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case has the p^-wave symmetry, and midgap states exist 
at the ends of the chain. They are described by the Ma- 
jorana fermions, which Kitaev proposed to use for non- 
volatile memory in quantum computing. He found that, 
when two such superconductors are brought in contact, 
the system is 47r-periodic in the phase difference between 
the superconductors. Our results are in agreement with 
his work. However, we formulate the problem as an ex- 
perimentally realistic Joscphson effect between known su- 
perconducting materials. 

For completeness, we also calculate the spectrum of 
Andreev bound states and the Josephson current between 
a singlet s-wave and a triplet p-wave superconductors, as 
well as between two chiral p-w&ve superconductors [28] . 
In agreement with previous literature [29-31], we find 
that a Josephson current is permitted between singlet 
and triplet superconductors, contrary to a common mis- 
conception that it is forbidden by the symmetry differ- 
ence. However, we do not find the fractional Josephson 
effect in these cases. 



II. THE BASICS 

The spin symmetry of the Cooper pairing is classified 
as either singlet (c CT (fc)<v(— fc)) oc e aa >A(k) = ia^,A(k) 
or triplet (c CT (fe)c^(-fe)) oc ia^(& ■ ra)A(fc) [32]. Here 
c CT (fc) is the annihilation operator of an electron with 
the spin a and momentum fc; e CTCT / is the antisymmet- 
ric metric tensor and <x are the Pauli matrices acting 
in the spin space; n is a unit vector characterizing po- 
larization of the triplet state. In this paper, we con- 
sider only the class of triplet superconductors where the 
spin-polarization vector n has a uniform, momentum- 
independent orientation. Everywhere in the paper, ex- 
cept in Sec. IIIF, we select the spin quantization axis z 
along the vector n. Then the Cooper pairing takes place 
between electrons with the opposite z-axis spin projec- 
tions a and u: (c (7 (k)c s (—k)} oc A a (k). Because the 
fermion operators c anticommute, the pairing potential 
has the symmetry A CT (fe) = ^fA s (k) = ±A CT (— k), where 
the upper and lower signs correspond to the singlet and 
triplet cases. 

We select the coordinate axis x perpendicular to the 
Josephson junction plane. We assume that the interface 
between the two superconductors is smooth enough, so 
that the electron momentum component k y , parallel to 
the junction plane, is a conserved good quantum number. 

Electron states in a superconductor are described by 
the Bogoliubov operators 7, which are related to the elec- 
tron operators c by the following equations [33] 



lnak y = J dx [U* naky (x) C aky (X) + V* naky (x) cj_^ (x)^) 

Ct*, ( x ) = J2 K-fc B fa) %ak y + v* n -- ky (x) jl- h ] , (2) 

n 

where k y — ~k y , and n is the quantum number of 



the Bogoliubov eigenstates. The two-components vec- 
tors VWfc,, (x) = [u naky (x) (x)] are the eigenstates 
of the Bogoliubov-de Gennes (BdG) equation with the 
eigenenergies E nrjky 



£ky(k x ) + U(x) A Gky (x, k x ) 

KkAx,k) -£k y (k)-u(x) 



1p n = E n 1p n , (3) 



where k x = —id x is the x component of the electron 
momentum operator, and U(x) is a potential. In Eq. (3) 
and below, we often omit the indices a and k y to shorten 
notation where it does not cause confusion. 



III. JUNCTIONS BETWEEN 
QUASI-ONE-DIMENSIONAL 
SUPERCONDUCTORS 

In this section, we consider junctions between two 
Q1D superconductors, such as organic superconductors 
(TMTSF)2X, with the chains along the x axis, as shown 
in Fig. 1(a). For a Q1D conductor, the electron energy 
dispersion in Eq. (3) can be written as e = h 2 k1/2m — 
2tbCOs(bk y ) — n, where m is an effective mass, \x is the 
chemical potential, b and tb are the distance and the tun- 
neling amplitude between the chains. The superconduct- 
ing pairing potentials in the s- and p x -w&ve cases have 
the forms 



Ao-fe B (2;, k x ) 



crAp, s-wave, 
Apk x /k F , p x -w&ve, 



(4) 



where HkF = \f2mfi is the Fermi momentum, and a is 
treated as + for j and — for J.. The index j3 = R,L labels 
the right (x > 0) and left (x < 0) sides of the junction, 
and Af) acquires a phase difference <p across the junction: 



At 



A e 



A, 



A . 



(5) 



The potential U{x) = UnS(x) in Eq. (3) represents the 
junction barrier located at x = 0. Integrating Eq. (3) 
over x from -0 to +0, we find the boundary conditions 
at x = 0: 



?pL =i>R, x 1pR - 

Z = 2mUo/h k F , 



■d x i> L = k F Zi)(0), (6) 
£ = 4/(Z 2 + 4), (7) 



where D is the transmission coefficient of the barrier. 



A. Andreev bound states 

A general solution of Eq. (3) is a superposition of the 
terms with the momenta close to akp, where the index 
a = ± labels the right- and left-moving electrons: 



O 0K X 



Uf3a+ 
Vf3a+ 



ikpx 



Br-, 



U0a- 



(8) 
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where (3 = =F for R and L. Eq. (8) describes a subgap 
state with an energy \E\ < A , which is localized at the 
junction and decays exponentially in x within the length 
The coefficients (up aa , vp aa ) in Eq. (8) are deter- 
mined by substituting the right- and left-moving terms 
separately into Eq. (3) for x ^ 0, where U(x) = 0. In the 
limit kp 3> k, we find 



(a) 



A 



Aoe" 



V/3o 



E + iaflfiKVp 

Ap aa 



hvp 



where v F — Tikp/m is the Fermi velocity, and 



(9) 



crAp, s-wave, 
aAp, p^-wave, 



(10) 



with A^ given by Eq. (5). The fcy-dependent Fermi mo- 
mentum hkp = hkp + 2tbCOs(bk y )/vp in Eq. (8) elimi- 
nates the dispersion in k y from the BdG equation. 

Substituting Eq. (8) into the boundary conditions (6), 
we obtain four linear homogeneous equations for the co- 
efficients Ap and Bp. These equations are compatible if 
the determinant of the corresponding 4x4 matrix is zero. 
This compatibility condition has the following form: 



- V +<T _U_ (T _)(u +cr+ V_ - + - V +rT+ U-^ 



(w +cr _u_ CT+ - u +(7 _it_ - + )(ti +0 - + t>_ (7 _ - v +a+ u- a -) 
= 1-D. (11) 

Using the variables r\ defined in Eq. (9), Eq. (11) can be 
written in a simpler form 



(V-<r- - »?+<t-)(»7-<t+ - V+v+) 



<7+ 



V+a-)(V-o- -V+a+) 



1 -D. 



(12) 



Substituting Eq. (9) into Eq. (12), we obtain an equation 
for the energies of the Andreev bound states. For a given 
<r, there are two subgap states with the energies E a = 
aE (<j>) labeled by the index a = ±, where 



E ( s) (<P) = -A yl -L>sin 2 (0/2), s-s junction,(13) 
Eo ] (4>) = -AoVDcoa{4>/2), p x -p x junction. (14) 



The energies (13) and (14) are plotted as functions 
of <f) in the left panels (b) and (c) of Fig. 1. Without 
barrier (D = 1), the spectra of the s-s and p x -p x junc- 
tions are the same and consist of two crossing curves 
E = =fA o cos0/2, shown by the thin lines in the left 
panel of Fig. 1(b). A nonzero barrier (D < 1) affects 
the energies of the Andreev bound states in the s-s and 
p x -p x junctions in different ways. In the s-s case, the 
two energy levels repel near <f> — tt and form two sepa- 
rated 27r-periodic branches shown by the thick lines in 
the left panel of Fig. 1(b). This is well known for the s-s 
junctions [33, 34]. In contrast, in the p x -p x case, the two 
energy levels continue to cross at <j> = 7r, and they detach 
from the continuum of states above +Ao and below — Ao 
at <f> — and 27r, as shown in the left panel of Fig. 1(c). 
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FIG. 1: (a) Butt-to-butt Josephson junction between two 
Q1D p^-wave superconductors, (b) The energies (left panel) 
and the currents (right panel) of the Andreev subgap states 
in the s-s junction as functions of the phase difference <fr for 
D = 1 (thin lines) and D = 0.9 (thick lines), (c) The same as 
(b) for the p x -p x junction at D — 0.2. 



The absence of energy levels repulsion at <f> = ir indicates 
that there is no matrix element between these levels in 
the p x -p x case, unlike in the s-s case. 

As shown in Sec. IV A, the 45°/45° junction between 
two d-wave superconductors is mathematically equiva- 
lent to the p x -p x junction. Eq. (14) was derived for the 
45°/45° junction in Refs. [24, 25, 35]. 



B. dc Josephson effect in thermodynamic 
equilibrium 

It is well known [33, 36] that the current carried by a 
quasiparticle state a is 



2edEa 
a ~ h 



(15) 



The two subgap states carry opposite currents, which 
are plotted vs. <p m the right panels (b) and (c) of Fig. 
1 for the s-s and p x -p x junctions. In thermodynamic 
equilibrium, the total current is determined by the Fermi 
occupation numbers f a of the states at a temperature T: 



T 2e ^ dE a 2e dE . . 

n = -\- ' 



h d<j) 



Eo 
2T 



(16) 



For the s-s junction, substituting Eq. (13) into Eq. (16), 
we recover the Ambegaokar-Baratoff formula [37] in the 
tunneling limit D <C 1 



eA /A N 
s D sin <b — — tanh — — 
Y 2H \2T y 



= Sin ^ tanh (^) 
(17) 
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FIG. 2: Critical currents of the s-s (dashed line) and p x -p x 
(solid line) Josephson junctions as functions of temperature 
for D = 0.3. 



and the Kulik-Omelyanchuk formula [38] in the transpar- 
ent limit D — > 1 



sin 



cos(0/2)\ 



2T 



(18) 



Taking into account that the total current is proportional 
to the number N of conducting channels in the junction 
(e.g. the number of chains), we have replaced the trans- 
mission coefficient D in Eq. (17) by the junction resis- 
tance R = h/2Ne 2 D in the normal state. 

Substituting Eq. (14) into Eq. (16), we find the Joseph- 
son current in the p x -p x junction in thermodynamic equi- 
librium: 



1 „ = V^Dsin ( — ) - tanh 



= sm 



t 
2 

ttA 



'a VD 



cos 



(0/2) \ 
2T J 



DeR 



tanh 



A Q VD cos(0/2) 
2T 



(19) 



The temperature dependences of the critical currents 
for the s-s and p x -p x junctions are shown in Fig. 2. They 
are obtained from Eqs. (17) and (19) assuming the BCS 
temperature dependence for Ao. In the vicinity of T c , 
I p and I s have the same behavior. With the decrease 
of temperature, I s quickly saturates to a constant value, 
because, for D < 1, E ( a s) w =fA (13), thus, for T < A , 
the upper subgap state is empty and the lower one is com- 
pletely filled. In contrast, I p rapidly increases with de- 
creasing temperature as 1 /T and saturates to a value en- 
hanced by the factor 2/\/~D relative to the Ambegaokar- 
Baratoff formula (13) at T = 0. This is a consequence 
of two effects. As Eqs. (17) and (19) show, I s oc D and 
I p oc >/D, thus I p S> I s in the tunneling limit Del. At 
the same time, the energy splitting between the two sub- 
gap states in the p x -p x junction is small compared to the 

gap: oc VdA < A . Thus, for y/DA < T < A , 
the two subgap states are almost equally populated, so 
the critical current has the 1/T temperature dependence 
analogous to the Curie spin susceptibility. 



Eq. (19) was derived analytically for the 45°/45° 
junction between two <i-wave superconductors in Refs. 
[23, 24], and a similar result was calculated numerically 
for the p x -p x junction in Ref. [18, 19]. Notice that Eq. 
(19) gives the Josephson current I p {<f>) that is a 27r- 
periodic functions of (f>, both for T — and T ^ 0. This 
is a consequence of the thermodynamic equilibrium as- 
sumption. At T = 0, this assumption implies that the 
subgap state with the lower energy is occupied, and the 
one with the higher energy is empty. As one can see in 
Fig. 1 , the lower energy is always a 27r-periodic functions 
of 4>. The assumption of thermodynamic equilibrium was 
explicitly made in Ref. [24] and was implicitly invoked 
in Refs. [18, 19, 23] by using the Matsubara diagram 
technique. In Ref. [39], temperature dependence of the 
Josephson critical current was measured in the YBCO 
ramp-edge junctions with different crystal angles and was 
found to be qualitatively consistent with the upper curve 
in Fig. 2. 



C. Dynamical fractional ac Josephson effect 

The calculations of the previous section apply in the 
static case, where a given phase difference <j> is maintained 
for an infinitely long time, so the occupation numbers of 
the subgap states have enough time to relax to thermo- 
dynamic equilibrium. Now let us consider the opposite, 
dynamical limit. Suppose a small voltage eV <C Ao is 
applied to the junction, so the phase difference acquires 
dependence on time t: <j)(t) — 2eVt/h. In this case, the 
state of the system is determined dynamically starting 
from the initial conditions. Let us consider the p x -p x 
junction at T = in the initial state <j> = 0, where the 
two subgap states (14) with the energies ±E are, cor- 
respondingly, occupied and empty. If <f>(t) changes suf- 
ficiently slowly (adiabatically), the occupation numbers 
of the subgap states do not change. In other words, the 
states shown by the solid and dotted lines in Fig. 1(c) 
remains, correspondingly, occupied and empty. The oc- 
cupied state (14) produces the current (15): 



Ip(t) 



sm 



DeAr 



sm 



eVt 



(20) 

The frequency of the ac current (20) is eV/%, a half of the 
conventional Josephson frequency 2eV/h. The fractional 
frequency can be traced to the fact that the energies Eq. 
(14) and the corresponding wave functions have the pe- 
riod 47r in <fi, rather than conventional 2-7T. Although at 
<t> = 2tt the spectrum in the left panel of Fig. 1(c) is the 
same as at <f> — 0, the occupation numbers are different: 
The lower state is empty and the upper state is occupied. 
Only at <j> = 4ir the occupation numbers are the same as 
at (f> = 0. 

The 47r periodicity is the consequence of the energy 
levels crossing at <f> — it. (In contrast, in the s-wave 
case, the levels repel at <f> = n in Fig. 1(b), thus the en- 
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ergy curves are 27r-periodic.) As discussed at the end of 
Sec. Ill A, there is no matrix element between the cross- 
ing energy levels at <fi — n - Thus, there are no transi- 
tions between them, so the occupation numbers of the 
solid and dotted curves in Fig. 1(c) are preserved. In 
order to show this more formally, we can write a gen- 
eral solution of the time-dependent BdG equation as a 
superposition of the two subgap states with the time- 
dependent <j>(i): ip(t) = X) C (i)^a [<£(*)]■ The matrix 
element of transitions between the states is proportional 
to = 4>{il}+\d <t> H\ij-)/{E- - E+). We found 

that it is zero in the p^-wave case, thus there are no tran- 
sitions, and the initial occupation numbers of the subgap 
states at <f> — are preserved dynamically. 

As one can see in Fig. 1(c), the system is not in the 
ground state when tt < <fi < 3tt, because the upper energy 
level is occupied and the lower one is empty. In princi- 
ple, the system might be able to relax to the ground state 
by emitting a phonon or a photon. At present time, we 
do not have an explicit estimate for such inelastic re- 
laxation time, but we expect that it is quite long. (The 
other papers [18, 19, 23, 24] that assume thermodynamic 
equilibrium for each value of the phase <f> do not have 
an estimate of the relaxation time either.) To observe 
the predicted ac Josephson effect with the fractional fre- 
quency eV/h, the period of Josephson oscillations should 
be set shorter than the inelastic relaxation time, but not 
too short, so that the time evolution of the BdG equation 
can be treated adiabatically. Controlled nonequilibrium 
population of the upper Andreev bound state was re- 
cently achieved experimentally in an s-wave Josephson 
junction in Ref. [40]. 

Eq. (20) can be generalized to the case where initially 
the two subgap states are populated thermally at = 0, 
and these occupation numbers are preserved by dynami- 
cal evolution 

m . | E ^M! /[E „ w= o)] ,2i, 



= sin 



h ^ d<j) 
'eVt\ ttA 



h ) \f~DeR 



««»h(^Y ,22, 



Notice that the periodicities of the dynamical equation 
(22) and the thermodynamic Eq. (19) are different. The 
latter equation assumes that the occupation numbers of 
the subgap states are in instantaneous thermal equilib- 
rium for each d>. 



D. Tunneling Hamiltonian approach 

In the infinite barrier limit D — > 0, the energies iE^ 
of the two subgap states (14) degenerate to zero, i.e. they 
become midgap states. The wave functions (8) simplify 
as follows: 



^ ±0 = 71 



(23) 



ipw = V 2k sin(fci? x) e K 
ipno = V%k sm(kpx) e~ 



9{-x), 



(24) 



i4>/2 

_l-i H 2 ) *(*)■ (25) 



Since at D — the Josephson junction consists of two 
semi-infinite uncoupled p^-wave superconductors, V'lo 
and ipRo are the wave functions of the surface midgap 
states [11] belonging to the left and right superconduc- 
tors. Let us examine the properties of the midgap states 
in more detail. 

If (u, v) is an eigenvector of Eq. (3) with an eigenvalue 
E n , then (— v*, u*) for s-wave and (v*,u*) for p-wave are 
the eigenvectors with the energy E n = —E n . It follows 
from these relations and Eq. (1) that ^ nsky = C^ ak 
with |C| = 1. Notice that in the s-wave case, because 
(u,v) and (— v*,u*) arc orthogonal for any u and v, the 
states n and n are always different. However, in the p- 
wave case, the vectors (u, v) and (v*,u*) may be propor- 
tional, in which case they describe the same state with 
E = 0. The states (24) and (25) indeed have this prop- 
erty: 



iu 



L0> 



-iu 



RO- 



(26) 



Substituting Eq. (26) into Eq. (1), we find the Bogoliubov 
operators of the left and right midgap states 



~ l lR0&k v 



(27) 



Operators (27) correspond to the Majorana fermions dis- 
cussed in Ref. [27]. In the presence of a midgap state, 
the sum over n in Eq. (2) should be understood as 
Sn>o"'"(V^) Sn=0' where we identify the second term 
as the projection Vc of the electron operator onto the 
midgap state. Using Eqs. (26), (27), and (2), we find 



'Pc aky (x) = u {x)% aky 



u o(*)7(L 



(28) 



Let us consider two semi-infinite p x -wave superconduc- 
tors on a ID lattice with the spacing I, one occupying 
x < I — —I and another x > I. They are coupled by the 
tunneling matrix element r between the sites I and /: 

Hr =r^Lk y (J)^k y Q) + 2L* S (0£W0]- (29) 

crky 

In the absence of coupling (t = 0) , the subgap wave func- 
tions of each superconductor are given by Eqs. (24) and 
(25). Using Eqs. (28), (26), (24), and (25), the tunneling 
Hamiltonian projected onto the basis of midgap states is 

VH T = t [u* L0 {l)u m {l) + ex.] [ll 0] lR0] + H.c.) 
= A Vd cos(0/2) (7l OT 7iioT + l R0 ^m), (30) 

where \J~D = 4r sin 2 kp l/Kvp is the transmission ampli- 
tude, and we omitted summation over the diagonal index 
k y . Notice that Eq. (30) is 47r-periodic in <j> [27]. 
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Hamiltonian (30) operates between the two degener- 
ate states of the system related by annihilation of the 
Bogoliubov quasiparticlc in the right midgap state and 
its creation in the left midgap state. In this basis, Hamil- 
tonian (30) can be written as a 2 x 2 matrix 



VH T = A o Vz?cos(0/2) 



1 

1 



(31) 



The eigenvectors of Hamiltonian (31) are (l,=Fl), i.e. the 
antisymmetric and symmetric combinations of the right 
and left midgap states given in Eq. (23). Their eigenen- 
ergies are E±(4>) = ^fA V~D cos(^/2), in agreement with 
Eq. (14). The tunneling current operator is obtained by 
differentiating Eqs. (30) or (31) with respect to <j>. Be- 
cause <f> appears only in the prefactor, the operator struc- 
tures of the current operator and the Hamiltonian are the 
same, so they are diagonal in the same basis. Thus, the 
energy eigenstates are simultaneously the eigenstates of 
the current operator with the eigenvalues 



/±=± ^^ Sm V2 



(32) 



in agreement with Eq. (20). The same basis (l,Tl) di- 
agonalizcs Hamiltonian (31) even when a voltage V is 
applied and the phase tfi is time-dependent. Then the 
initially populated eigenstate with the lower energy pro- 
duces the current I p = \/~D(eAo/h) sui(eVt/%) with the 
fractional Josephson frequency eV/Ti, in agreement with 
Eq. (20). 



E. Josephson current carried by single electrons, 
rather than Cooper pairs 

In the tunneling limit, the transmission coefficient D 
is proportional to the square of the electron tunneling 
amplitude r: D oc r 2 . Eqs. (20) and (32) show that the 
Josephson current in the p x ~Px junction is proportional 
to the first power of the electron tunneling amplitude t. 
This is in contrast to the s-s junction, where the Joseph- 
son current (17) is proportional to t 2 . This difference 
results in the big ratio I p /I s = 2/%/!? between the crit- 
ical currents at T — in the p x - and s-wave cases, as 
shown in Fig. 2 and discussed in Sec. IIIB. The reason 
for the different powers of r is the following. In the p x - 
wave case, the transfer of just one electron between the 
degenerate left and right midgap states is a real (non- 
virtual) process. Thus, the eigenenergies are determined 
from the secular equation (31) already in the first order 
of r. In the s-wave case, there are no midgap states, so 
the transferred electron is taken from below the gap and 
placed above the gap, at the energy cost 2Ao. Thus, the 
transfer of a single electron is a virtual (not real) process. 
It must be followed by the transfer of another electron, so 
that the pair of electrons is absorbed into the condensate. 
This gives the current proportional to r 2 . 



This picture implies that the Josephson supercurrcnt 
across the interface is carried by single electrons in the p x - 
p x junction and by Cooper pairs in the s-s junction. Be- 
cause the single-electron charge e is a half of the Cooper- 
pair charge 2e, the frequency of the ac Josephson effect 
in the p x -p x junction is eV/h, a half of the conventional 
Josephson frequency 2eV/h for the s-s junction. These 
conclusions also apply to a junction between two cuprate 
ci-wave superconductors in such orientation that both 
sides of the junction have surface midgap states, e.g. to 
the 45°/45° junction (see Sec. IV A). 

In both the p x -p x and s-s junctions, electrons trans- 
ferred across the interface are taken away into the bulk by 
the supercurrent of Cooper pairs. In the case of the p x -p x 
junction, a single transferred electron occupies a midgap 
state until another electron gets transferred. Then the 
pair of electrons becomes absorbed into the bulk conden- 
sate, the midgap state returns to the original configura- 
tion, and the cycle repeats. In the case of the s-s junc- 
tion, two electrons are simultaneously transferred across 
the interface and become absorbed into the condensate. 
Clearly, electric charge is transferred across the interface 
by single electrons at the rate proportional to r in the 
first case and by Cooper pairs at the rate proportional 
to t 2 in the second case, but the bulk supercurrcnt is 
carried by the Cooper pairs in both cases. 



F. Josephson effect between triplet 
superconductors with nonparallel n-vectors 

In this section, we consider the Josephson effect be- 
tween two p x -wave superconductors with nonparallel 
spin-polarization vectors n forming an angle 9. This 
problem was studied in Ref. [19] using a tunneling Hamil- 
tonian approach. Here we analyze the problem using 
the BdG formulation. There are experimental indica- 
tions that the spin-polarization vector n is parallel to 
the crystal axis c in the (TMTSF) 2 X compounds [3, 5]. 
Then the considered junction can be realized in the ge- 
ometry shown in Fig. 1(a) where the c axes of the two 
superconductors are rotated relative to each other by the 
angle 9 around the common a axis along the chains. 

Let us select the spin quantization axis z perpendicular 
to both vectors n, and the x axis in the spin space parallel 
to the vector n of the left superconductor. Then the 
vector n of the right superconductor lies in the (x,y) 
plane at the angle 9 to the x axes: n = (cos 9, sin 9,0). 
In this representation, the superconducting pairing takes 
place between electrons with parallel spins: 



(c a (k)c a ,(-k)) cx ia {y) {a {x) n x + a {v) n y )A{k) 

A(fe). (33) 



e lV 
-e" 



Then, the Josephson effect can be considered separately 
for the spin up and down sectors having the phase dif- 
ferences (j) ± 9, correspondingly. Using Eq. (14) for the 
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Px-Px junction, we find the energies of the Andreev bound 
states for each spin sector 



E a , 



E a ,i 



-aA \/D cos 



-aA vD cos 



(34) 



(35) 



The total Josephson current is obtained by adding the 
currents carried by the two spin sectors [41]. For simplic- 
ity, below we consider only the case of zero temperature. 
In the dynamical limit, assuming that the states (34) and 
(35) with a = + are occupied initially and the occupation 
numbers are preserved dynamically and using Eq. (20), 
we find a 47r-periodic current: 



I(t) = 



VDeAp 

h 



sm 



sm I — — I cos 



(36) 



In the static thermodynamic limit, using Eq. (19) at T = 
0, we find the dc Josephson current: 



DeA a 
2h 



sm 



sgn 



cos 



sm 



sgn 



cos 



(37) 



For completeness, let us also consider the Josephson 
effect between two p y -w&ve or two p z -wave superconduc- 
tors, where the y and z axes are parallel to the junction 
plane. In these junctions, midgap states are absent in the 
D — > limit, thus the current-phase relation is conven- 
tional I = I c sin 4>. For nonparallel vectors n, the total 
Josephson current is the sum of the spin up and down 
sectors: 



y [sin(0 + 
I r cos 6 sin i 



9) + sin(0 - 9)} 
b = I c {n L ■ n R )sin( 



(38) 



Eq. (38) is consistent with Ref. [19] . In the case where the 
two vectors n are perpendicular (9 = ir/2), the Joseph- 
son current (38) for the superconductors without midgap 
states vanishes, but, according to Eqs. (36) and (37), it 
is not zero if the midgap states are present. 



G. s-p x junction between singlet and triplet 
superconductors 



In this section, we consider a junction between a singlet 
s-wave and a triplet p x -wave superconductors. The junc- 
tion geometry is the same as in Fig. 1(a), where one of 
the superconductors is taken to be a conventional s-wave 
superconductor and another one a Q1D triplet p^-wave 
superconductor, such as (TMTSF) 2 X. 



(a) 1 
0.5 

< 

LU 

-0.5 





FIG. 3: The subgap energy levels and the Josephson currents 
in the s-p x junction. Here D = 0.8 and Io = eAo/2%. (a) 
The energies (39) of the Andreev bound states. The solid and 
dashed lines correspond to a =| and a =|. (b) The Josephson 
current in the static thermodynamic limit, where the states 
with E < are occupied, (c) The Josephson current in the 
dynamical limit, where the central branch with a = j and the 
lower branches touching — Ao are occupied. 



We choose the spin quantization axis z along the po- 
larization vector n of the triplet superconductor, so the 
spin projection a on the z axis is a good quantum num- 
ber. In both triplet and singlet superconductors, the 
Cooper pairing takes place between electrons with oppo- 
site spins. However, the pairing potential has the same 
sign for a and a in the triplet superconductor and the 
opposite signs in the singlet superconductor. Thus, the 
phase difference across the Josephson junction is <f> for 
quasiparticles with a =\ and cj) + n for a —\. The ener- 
gies of the Andreev bound states can be found for each a 
from Eq. (12) together with Eq. (9), where we should use 
the upper line of Eq. (10) for the left superconductor and 
the lower line for the right superconductor. To simplify 
calculations, we consider the case where the magnitudes 
of the gaps are equal for the s- and p^-wave superconduc- 
tors: | At | = |Afl| = A . The energies of the Andreev 
bound states are 



E a ,a = a- sgn (sin 0)A O 1 



1 + ay/1 - D 2 sin 2 



(39) 



For each value of the spin index a = ±, Eq. (39) gives 
two Andreev states labeled by the index a = ±. In the 
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tunneling limit D <C 1, we have 



E 



E. 



a sgn(sin < 



- aA D sin( 



1 - 



£ 2 sin 2 



cj>), (40) 
(41) 



The energies (39) are plotted in Fig. 3(a) vs. <fi by the 
solid lines for a —\ and by the dashed lines for a ={. Wc 
observe that the branches (40) with a = + touch the gap 
boundaries ±A at = and it, whereas the branches 
(41) with a — — stay in the center of the gap. 

In the limit D C 1, the central branches with a = — 
dominate the energy dependence on <f>, and energy min- 
ima are achieved at <j> — 7T /2 or 37r/2. Notice that if 
the system selects the energy minimum at <f> = tt/2, then 
the spin down states, shown by the dashed lines in Fig. 
3(a), are populated, and the spin up states are empty, so 
the junction accumulates the spin —h/2 per conducting 
channel [41]. If the system selects the energy minimum 
at <j> = 3ir/2, then the junction accumulates the spin h/2 
per conducting channel. 

In the limit D <C 1, we can neglect the energies (40) 
and obtain the Josephson current by differentiating the 
energies (41) with respect to (f> using Eq. (15) [41]. In the 
dynamical limit, the occupation numbers of the Andreev 
states (41) are preserved, and the Josephson current has 
the 27r-periodicity, as shown in Fig. 3(c): 



J(t) 



2h 



A Dcos(/)(t). 



(42) 



In the static thermodynamic limit, the system occupies 
the branch of the minimal energy for each <j>, and the 
Josephson current is 7r-periodic, as shown in Fig. 3(b): 



I(4>) ~ -^A L>sgn(sin< 



COS ( 



(43) 



The thermodynamic assumption implies that the spin ac- 
cumulation at the s-p x junction changes sign when the 
phase <p crosses w. 

Now let us consider the circuit shown in Fig. 4, where a 
p x -wave superconductor has Josephson junctions at both 
ends with an s-wave superconductor closed in a loop. Be- 
cause the sign of the p^-wave pairing potential is opposite 
for the +kp and — kp sheets of the Fermi surface, the 
two junctions have the relative phase shift ir. Naively, 
one might expect a spontaneous current in this circuit 
by analogy with the corner SQUID in the cuprates [1]. 
However, the system shown in Fig. 4 can accommodate 
the phase shift tt by selecting the energy minimum at 
(f> = n/2 for one junction and the energy minimum at 
<f> = 3tt/2 for another junction. Then, no current cir- 
culates in the loop. However, one junction accumulates 
spins up and another junction spins down, which might 
be possible to detect experimentally. 

The results of the this section clearly show that a 
Josephson current is possible between singlet and triplet 
superconductors, in agreement with the earlier findings 




s-wave 



FIG. 4: A Q1D p^-wave superconductor closed in a loop by 
an s-wave superconductor. No current is circulating in the 
loop in equilibrium. However, there is accumulation of spins 
up in one s-p x junction and spins down in another junction. 
The sketch at the top illustrates the Fermi surface of a Q ID 
metal with the opposite signs of the superconducting p^-wave 
pairing potential on the two sheets of the Fermi surface. 



by Yip [29]. Recently, the Josephson current was calcu- 
lated for the s-p x junction in Ref. [31], but spin accumu- 
lation at the junction was not recognized in this paper. 
The s-px junction considered in this section is mathemat- 
ically equivalent to the 0°/45° d-d junction and the 45° 
junction between an s-wave and a d-w&ve superconduc- 
tors (see Sec. IV A). Eq. (39) was obtained for that case 
in Rcfs. [23-25]. However, there is no spin accumulation 
in junctions between singlet s- and <i-wave superconduc- 
tors, unlike in the s-p x junction. 



IV. JUNCTIONS BETWEEN 
QUASI-TWO-DIMENSIONAL 
SUPERCONDUCTORS 

In this section, we study junctions between quasi-two- 
dimensional (Q2D) superconductors such as nonchiral d- 
wave cuprates and chiral p x ± ip y -wave ruthenates. For 
simplicity, we use an isotropic electron energy dispersion 
law e = fi 2 (kx + ky)/2m — fi in the (x,y) plane. As 
before, we select the coordinate x perpendicular to the 
junction line and assume that the electron momentum 
component k y parallel to the junction line is a conserved 
good quantum number. Then, the 2D problem separates 
into a set of ID solutions (8) in the x direction labeled 
by the index k y . The Fermi momentum kp and velocity 

vp are replaced by their x-components kp x = 



kp 



k y 



and vp x — hkp x /m. The transmission coefficient of the 
barrier (7) becomes fc y -dependent 



Z{k v ) = Z Q 



kr 



lh 2 — £-2 
Y F S 



D{ky) = 



(44) 



where Z = ^2mU /fikp. The total Josephson current 
is given by a sum over all occupied subgap states labeled 
by k y . 
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FIG. 5: Schematic drawing of the 45°/45° junction (panel a) 
and 0°/0° junction (panel b) between two d-wave supercon- 
ductors. The thick line represents the junctions line. The 
circles illustrate the Fermi surfaces, where positive and neg- 
ative pairing potentials A are shown by the solid and dotted 
lines. The points A, B, C, and D in the momentum space are 
connected by transmission and reflection from the barrier. 



A. Josephson junctions between d-wave 
superconductors 

For the cuprates, let us consider a junction parallel to 
the [1,1] crystal direction in the (a, 6) plane and select 
the x axis along the diagonal [1, 1], as shown in Fig. 5(a). 
In these coordinates, the d-wave pairing potential is 

Arf 9 {x, k x ) = a2A k y k x /kp, (45) 

where the same notation as in Eq. (4) is used. Direct 
comparison of Eqs. (45) and (4) demonstrates that the 
d-wave superconductor with the 45° junction maps to 
the p x -wave superconductor by the substitution A — > 
a2Aok y /kp. Thus, the results obtained in Sec. Ill for the 
Px-Px junction apply to the 45°/45° junction between two 
d-wave superconductors with the appropriate integration 
over k y . The energies of the subgap Andreev states are 
given by Eq. (14) with the /Cy-dependent parameters Ao 
and D, and the energies and the wave functions are 47r- 
periodic functions of <j>. Thus, the ac Josephson current 
in the dynamical limit is 47r-periodic and has the frac- 
tional frequency eV/h, as in Eqs. (20), (22), and (32). 
The energies (14) of the subgap states [24, 35] and the 
dc Josephson current (19) in the thermodynamic limit 
[23, 24] were calculated for the 45°/45° d-d junction be- 
fore. However, these papers did not recognize the frac- 
tional, 47r-periodic character of the Josephson effect in 
the dynamical limit. 

On the other hand, if the junction is parallel to 
the [0,1] crystal direction, as shown in Fig. 5(b), then 
A a k y (x, k x ) = aAp (k^ — ky)/k 2 F . This pairing potential 

is an even function of k x , thus it is analogous to the s- 
wave pairing potential in Eq. (4). Then, the 0°/0° junc- 
tion between two d-wave superconductors is analogous 
to the s-s junction. It should exhibit the conventional 
27T-periodic Josephson effect with the frequency 2eV/K. 

For a generic orientation of the junction line, the d- 
wave pairing potential is p^-like for some momenta k y 



and s-like for other k y . Thus, the total Josephson current 
is a sum of the unconventional and conventional terms: 

/ = C sin(0/2) + C x sin(0) + . . . , (46) 

where Co and Ci are some coefficients. We expect that 
both terms in Eq. (46) are present for any real junction 
between d-wave superconductors because of imperfec- 
tions in junction orientation. However, the ratio Cq/Ci 
should be maximal for the junction shown in Fig. 5(a) 
and minimal for the junction shown in Fig. 5(b). In gen- 
eral, whenever the superconductors on both sides of the 
junction have surface midgap states, we expect to observe 
the 47r-periodic fractional ac Josephson effect. In princi- 
ple, the effect may be spoiled by the gapless quasiparti- 
cles that exist near the gap nodes in a d-wave supercon- 
ductor. However, they would affect only a small portion 
of the Fermi surface near the nodes, and the 47r-periodic 
Josephson effect should survive on the other parts of the 
Fermi surface, where the gap is big. 

The 45°/45° junction shown in Fig. 5(a) should not 
be confused with the 0°/45° d-d junction [42] or the 45° 
s-d junction [43, 44], much discussed in literature. None 
of the papers [42-44] treated the problem correctly, be- 
cause they did not take into account the Andreev bound 
states in the junction properly. The correct energy spec- 
trum of the Andreev bound states was obtained in Rcfs. 
[24, 25, 35]. In the 0°/45° d-d and 45° s-d junctions, only 
one superconductor has midgap states, thus these junc- 
tions are mathematically analogous to the s-p x junction 
considered in Sec. Ill G. The spectrum of the Andreev 
bound states is given by Eq. (39) without the factor er, 
because both superconductors are singlet. The energy 
levels are plotted vs. <\> in Fig. 3(a), where the solid and 
dashed lines represent not spin, but positive and nega- 
tive momenta k y . The junction has two energy minima 
at (f> = 7r/2 or 37r/2, where the states with only negative 
or positive momenta k y are filled, thus there are persis- 
tent currents along the junction line [45, 46]. (On the 
other hand, there is no spin accumulation, unlike in the 
s-Px junction discussed in Sec. IIIG.) In the thermo- 
dynamic limit, the current-phase relation shown in Fig. 
3(b) is 7r-periodic; however, it requires reversing the cur- 
rents along the junction line when (f> passes through 
or 7r. In the dynamical limit, the current-phase relation 
shown in Fig. 3(c) is 27r-periodic. The first two harmonics 
I = C\ sin(0) + Ci sin(2</>) have been recently observed 
experimentally in the 0°/45° d-d junction [47]. 

B. Josephson junctions with chiral superconductors 

In this section, we study junctions between the chiral 
p x ± ipy-w&ve superconductors Sr 2 Ru04, where the pair- 
ing potential is assumed to be A(k) = Ao(k x ± ik y )/kp 
[6], and the two signs correspond to opposite chiralities. 
We assume a uniform orientation of the spin-polarization 
vector n across the junction. This problem was inves- 
tigated in Ref. [28] using the Eilenberger equation for 
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Green's functions. It was found that the chiral subgap 
states at the junctions enhance the low-temperature crit- 
ical Josephson current in symmetric junctions. Here we 
use the BdG equation to obtain the spectrum of the An- 
dreev bound states. As before, we assume that the mo- 
mentum component k y parallel to the junction is con- 
served. Thus, the problem separates into a set of ID 
solutions in the x direction perpendicular to the junction 
plane, and we can use the method of Sec. Ill A. 

First we consider a junction between two superconduc- 
tors with opposite chiralities, as illustrated in the first 
column of Fig. 6(a). In this case, = Ao(t +ik y )/kF 
and Ar — e l ^Ao(k x — ik y )/kF- When the barrier is 
not transparent (D = 0), each superconductor has chiral 
Andreev surface states with the same energy dispersion 
E(k y ) = kyAo/kp [13]. The electron tunneling ampli- 
tude r oc y/D produces a matrix element mixing the two 
states in the first-order degenerate perturbation theory. 
Thus, the two energy spectra repel with the splitting pro- 
portional to y/T). From Eq. (12), we find the following 
subgap energies: 




k% -k y VI) cos | J . 

(47) 

The energy levels splitting oscillates with the period 47T 
as a function of 6: 8E 



4>: SE = (A /k F )^k F - k%VD cos <j>/2 



The splitting depends on k y through the square-root pref- 
actor and through the dependence of D on k y in Eq. (44) , 
and vanishes at k y = ±kp. The energy dispersion (47) 
is plotted vs. k y in the second and third columns of Fig. 
6(a) for several values of <j). The spectrum of excitations 
is gapless because of the chiral dispersion in k y . Thus, 
it is reasonable to assume that the occupation numbers 
of the subgap states are in instantaneous thermodynamic 
equilibrium for any phase (f>. Then, the Josephson current 
is a 27r-periodic function of <fi, as illustrated at zero tem- 
perature in the fourth column of Fig. 6(a), even though 
the energy levels (47) are 47r-periodic functions of <f>. 

Now let us consider the case of two superconductors 
with opposite chiralities, as illustrated in the first column 
of Fig. 6(b). When the two superconductors are discon- 
nected (D = 0), their chiral Andreev surface states have 
opposite dispersions E = ±k y Ao/kp, thus they are non- 
degenerate. The electron tunneling amplitude r cx \[T) 
repels the energy levels around the intersection point 
k y = 0. From Eq. (12), we find the following subgap 
energies: 



E = ±-^\j{l-D)kl+ Dk% cos 2 



(48) 



The energy dispersion (48) is plotted vs. k y in the sec- 
ond and third columns of Fig. 6(b) for several values of 
<j>. The energy splitting around k y = is a 27r-periodic 
function of <j) and vanishes at 4> = ir. The Josephson cur- 
rent is a 27r-periodic function of <f>, as illustrated at zero 
temperature in the fourth column of Fig. 6(b). 



Now let us consider a junction between an s-wave and a 
Px + ipy-w&ve superconductors shown in the first column 
of Fig. 6(c). The Josephson current was calculated in this 
case in Ref. [31] using the method of Green's functions. 
However, the energies of the Andreev bound states were 
not written explicitly. The subgap states in this junction 
are obtained by solving Eq. (12) in the manner similar to 
the ID s-p x junction. For simplicity, we assume that the 
magnitudes of the pairing potentials in both supercon- 
ductors are the same: |Al| = |A#| = A . The square of 
the subgap energies is given by the following expression 

<a = ^-{l + Rkl-aDkySi^ (49) 
+ ay/l-kly/l- {Rky -aDsm(j))A , 



where R = 1 — D is the reflection coefficient, and k,„ 



ky/kF- The signs of the energies are 



sgni? a:(T = sgn I — Rk y + aD sin < 



(50) 



+ ak. 



\ 



{Rky - aD sin <^ 2 
\-~kl 



For a given a, there are two branches of energies labeled 
by the index a = ±. The energy dispersions E±^(k y ) 
are shown in the second and third columns of Fig. 6(c) 
for several phases 4>. In the limit of impenetrable barrier 
D — > 0, the energy branch with a = + approaches to the 
gap edges \E + \ — ► A , whereas the branch with a = — 
approaches to the energy dispersion — > —k y Ao/kF of 
the chiral surface states in the p x + zp^-wave supercon- 
ductor [13]. 

The energy E a ^ a (k y ) for a given a is a 27r-periodic func- 
tion of (f>. The energy E a . s {k y ) is obtained from E a ^ a {k y ) 
by the shift 4> ^ 4> + n , as discussed in Sec. Ill G. Thus, 
the Josephson current (16) in the static thermodynamic 
limit, obtained by summation over a and k yi is a tt- 
pcriodic function of cfi, as shown in the fourth column 
of Fig. 6(c), in agreement with Ref. [31]. Similarly to the 
s-p x junction considered in Sec. IIIG, the s-(p x + ip y ) 
junction has two equal energy minima [31] at <j) = t/2 
and 37r/2 accompanied by accumulation of the down spin 
for = 7r/2 and the up spin for <fi = 3n/2. 



V. EXPERIMENTAL OBSERVATION OF THE 
FRACTIONAL AC JOSEPHSON EFFECT 

Conceptually, the setup for experimental observation 
of the fractional ac Josephson effect is straightforward. 
One should apply a dc voltage V to the junction and 
measure frequency spectrum of microwave radiation from 
the junction, expecting to detect a peak at the fractional 
frequency eV/h. Higher harmonics, such as 2eV/h, may 
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also be present because of Eq. (46) and circuit nonlincar- 
ities, but an observation of the 1/2 subharmonic of the 
conventional Josephson frequency 2eV/h would be the 
signature of the effect. 

Josephson radiation at the conventional frequency 
2eV/h was first observed experimentally almost 40 years 
ago in Kharkov [48, 49], followed by further work [50, 51]. 
In Ref. [49], the spectrum of microwave radiation from 
tin junctions was measured, and a sharp peak at the fre- 
quency 2eV/h was found. Without any attempt to match 
impedances of the junction and waveguide, Dmitrenko 
and Yanson [49] found the signal several hundred times 
stronger than the noise and the ratio of linewidth to the 
Josephson frequency less than 10~ 3 . More recently, a 
peak of Josephson radiation was observed in Ref. [52] in 
indium junctions at the frequency 9 GHz with the width 
36 MHz. In Ref. [53], a peak of Josephson radiation 
was observed around 11 GHz with the width 50 MHz 
in Bi 2 Sr 2 CaCu208 single crystals with the current along 
the c axis perpendicular to the layers. 

To observe the fractional ac Josephson effect predicted 
in this paper, it is necessary to perform the same exper- 
iment with the 45°/45° cuprate junctions shown in Fig. 
5(a). For control purposes, it is also desirable to measure 
frequency spectrum for the 0°/0° junction shown in Fig. 
5(b), where a peak at the frequency eV/ft, should be min- 
imal. It should be absent completely in a conventional 
s-s junction, unless the junction enters a chaotic regime 
with period doubling [54, 55]. The high-T c junctions of 
the required geometry can be manufactured using the 
step-edge technique. Bicrystal junctions are not appro- 
priate, because the crystal axes a and b of the two su- 
perconductors are rotated relative to each other in such 
junctions. As shown in Fig. 5(a), we need the junction 
where the crystal axes of the two superconductors have 
the same orientation. Unfortunately, attempts to man- 
ufacture Josephson junctions from the Q1D organic su- 
perconductors (TMTSF) 2 X failed thus far. 

The most common way of studying the ac Josephson ef- 
fect is observation of the Shapiro steps [56]. In this setup, 
the Josephson junction is irradiated by microwaves with 
the frequency u>, and steps in dc current are detected 
at the dc voltages V n — nfiuj/2e. Unfortunately, this 
method is not very useful to study the effect that we pre- 
dict. Indeed, our results are effectively obtained by the 
substitution 2e — > e. Thus, we expect to see the Shapiro 
steps at the voltages V m = mhej/e = 2mhuj/2e, i.e. we 
expect to see only even Shapiro steps. However, when 
both terms are present in Eq. (46), they produce both 
even and odd Shapiro steps, so it would be difficult to dif- 
ferentiate the novel effect from the conventional Shapiro 
effect. Notice also that the so-called fractional Shapiro 
steps observed at the voltage V\ii — fiuj/Ae correspond- 
ing to n = 1/2 have nothing to do with the effect that 
we propose. They originate from the higher harmonics 
in the current-phase relation / oc sin(2</>). The fractional 
Shapiro steps have been observed in cuprates [57-59] , but 



also in conventional s-wave superconductors [60]. An- 
other method of measuring the current-phase relation in 
cuprates was employed in Ref. [61], but connection with 
our theoretical results is not clear at the moment. 



VI. CONCLUSIONS 

In this paper, we study suitably oriented p x -p x or d-d 
Josephson junctions, where the superconductors on both 
sides of the junction originally have the surface Andrccv 
midgap states. In such junctions, the Josephson cur- 
rent /, carried by the hybridized subgap Andrecv bound 
states, is a 47r-periodic function of the phase difference 
4>: I oc sin(</>/2), in agreement with Ref. [27]. Thus, the 
ac Josephson current should exhibit the fractional fre- 
quency eV/h, a half of the conventional Josephson fre- 
quency 2eV/h. In the tunneling limit, the Josephson 
current is proportional to the first power of the electron 
tunneling amplitude, not the square as in the conven- 
tional case [23-25]. Thus, the Josephson current in the 
considered case is carried across the interface by single 
electrons with charge e, rather than by Copper pairs with 
charge 2e. The fractional ac Josephson effect can be ob- 
served experimentally by measuring frequency spectrum 
of microwave radiation from the junction and detecting 
a peak at eV/h. 

In p x -p x junctions with nonparallel orientation of the 
spin-polarization vectors n, the Josephson current de- 
pends on the relative angle between the vectors n [19]. 
The Josephson current is permitted between singlet and 
triplet superconductors, but, in the static thermody- 
namic limit, the current-phase relation is 7r-pcriodic [29]. 
The s-p junction has two equal minima in energy at 
4> = tt/2 and 37r/2 [31], characterized by accumulation 
of the up or down spins (oriented relative to the vector 
n) in the junction. In Josephson junctions between chi- 
ral p x ± zpy-wave superconductors, the Andreev bound 
states are also chiral. In the static thermodynamic limit, 
the current-phase relation has the period of 27r in the 
chiral p-p junctions [28] and the period of it in the chiral 
s-p junctions [31]. 
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Note added in proof. The fractional Josephson ef- 
fect discussed in our paper is similar to the fractional 
quantum Hall effect [62]. Both involve existence of sev- 
eral equivalent ground states, whose energy levels cross: 
Compare Fig. lc of our paper and Fig. 2a of Ref. [62]. 
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FIG. 6: The subgap energy spectra and the Josephson cur- 
rents in the junctions with chiral p x ± ip H -wave superconduc- 
tors, calculated for Zq = 1. (a) Junction between two p-wave 
superconductors with opposite chiralities. (b) Junction be- 
tween two p-wave superconductors with the same chirality. 
(c) Junction between s-wave and chiral p-wave superconduc- 
tors. In the second and third columns of row (a), the solid 
and dashed lines show the a = ± branches of Eq. (47) for 
different values of <j). In rows (b) and (c), the energies (48) 
and (49) are shown by the solid and dashed lines for = 
and 7r/2 in the second column and for <f> = n and 37r/2 in the 
third column. The energies of quasiparticles with a =\ are 
the same as with a =f for rows (a) and (b), and they can 
be obtained by the shift <f> — » <f> + n for row (c). The fourth 
column shows the Josephson current in the static thermody- 
namic limit, normalized by Iq — eAoL y kF /2nh, where L y is 
the length of the junction. 



